%uitval transfo

%per bus kijken wat de spanning doet bij een uitval van een lijn.
%nu nog zonder vermogenvariatie

clear all
tic
IEEE_14_bus_original;
einde_interval = size(Line.con,1);%# bussen bekijken
teller = 0;
tijd = zeros(200,einde_interval);
for k = 1:size(Line.con,1)
    if Line.con(k,7)~=0
        start = k;
        break;
    end
end

aantal_trafos = (size(Line.con,1)-k+1);

alle_spanningen = zeros(200*aantal_trafos,14);

for k=start:einde_interval;
    
    filename = sprintf('trafofout%03d',k);
    initpsat;
    Settings.plottype = 2;
    Settings.freq = 60;
    Settings.absvalues= 'on';
    runpsat(filename,'pf');
    Settings.fixt = 1;
    Settings.tstep = 1;%0.05; %TRAAAAAAAAAAAAAG
    
    Settings.freq = 60;
    Settings.absvalues= 'on';
    Settings.plottype = 2;
    runpsat('td');
    
    %begint vanaf index 0... (STOM!!!)
   v = genvarname('voltages_fout_bus', who);
   eval([v ' = Varout.vars']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   lengte = size(Varout.vars,1);
   
   %spanning per tijdsstap per bus
   alle_spanningen(1+200*teller:200*teller+lengte,:) = Varout.vars(:,1:14);
   
   t = genvarname('tijd_fout_bus', who);
   eval([t ' = Varout.t']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   
   %tijd per stap per lijn
   tijd(1:lengte,k) = Varout.t;
   closepsat;
   teller = teller+1;
   
   save trafofout tijd alle_spanningen aantal_trafos;
end


